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Abstract 

A new approach l.o find homogeneous models for beam-like repeated 
flexible structures is proposed which conceptually involves two steps. 

Step one: Approximation of 3-D lion-homogeneous model by a 
1-D periodic beam model. The structure is modeled as a 3-D non- 
homogeneous continuum. The displacement field is approximated by 
lay lor series expansion. Then, the cross sectional mass and stiffness 
matrices are obtained by energy equivalence using their additive prop- 
erties. Due to the repeated nature of the flexible bodies, the mass and 
stiffness matrices are also periodic. This procedure is systematic and 
requires less dynamics detail. 

Step two: Homogenization from i-D periodic beam model to 1-D 
homogeneous beam model. The periodic beam model is homogenized 
into an equivalent homogeneous beam model using the additive prop- 
erty of compliance along the generic axis. The major departure from 
previous approaches in literature is using compliance instead of stiff- 
ness in homogenization. An obvious justification is that the stifFness 
is additive at each cross section but not along the generic axis. The ho- 
mogenized model preserves many properties of the original periodic 
model. 
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1 Introduction 

As the number of repeated cells in a truss structure increases, the 3-D model can 
be approximated better and better by an equivalent 1-D model. The repeated 
structure then can be modeled as a homogeneous anisotropic continuum beam. 
The parameters of the continuum beam are functions of the element properties of 
the truss structure. 

Finding the 1-D homogeneous anisotropic beam model from the reference model, 
the 3-D non-hoinogeneous anisotropic model, of the truss structure may be releired 
to as a homogenization process, of which there are many examples. 

The approach presented here follows that of Noor’s [1, 2] and Lions s [3j. See [3] 
for mathematical details of the homogenization process, where some results are taken 
by our paper for granted. The Noor’s method is a direct averaging method, which 
justifies equivalence in the sense of equal kinetic energy and potential energy under 
the condition of equal nodal displacements and velocities. It imposes a kinematic 
assumption on the displacement field, then averages the stiffness and mass matrices 
(by FEM) over a repeated cell. Although the stiffness matrix is additive at each 
cross-section, it. is not along the generic axis. Thus, this method always gives higher 
stiffness than it should be. This shortcoming will be overcome by our approach. 

Our approach consists of two steps, as illustrated in Figure (2). 1 he first 

step deals with the approximation of 3-D noil-homogeneous model of a repeated 
structure by the 1-D periodic beam model. The second step then homogenized 
it to a 1-D homogeneous beam model. The 3-D non-homogeneous model is a 
collection of the Eulerian Equation of Motion of each element of the structure, 
and is referred to as a reference model for the successive approximation. By applying 
the Taylor series expansion and energy equivalence, a 1-D periodic beam model is found 
systematically. Solid beam is used to clarify the basic idea, then an extension from 
solid beam to non-solid structure(eg. lattice structure) is presented in section (5). 



Figure 1: Anisotropic Deam 



Consider n structure constructed by linear elastic anisotropic materials as shown 
in figure (1). The material coordinates attached have x-axis as the genetic axis 
along the centroids, and y-z as the principal axes of the area inertia of cross 
sections. This choice of reference will be adopted throughout the paper. 

Let the bounded open set ft c r R denote the space occupied by the structure 
and T the boundary of ft. Let U, V, and W be the displacements in x, y and z 
direction, respectively, measured w.r.t. the natural state (undeformed position), 
observed from inertial coordinates, and represented in the material coordinates. 
The equntion of motion [4, 5] is 


pu = 

ftxx,x 4* &yx,y 4 * &zx,z 4 f r 



pv = 

ftxy,x 4 CTyj ,,y 4 Vzy,z 4 fy 
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with well-posed initial and boundary conditions to render existence of unique solu- 
tion, where 

P - P(x,y,z) € 

and f(x,y,z) is the external body force. The constitutive law is 
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where cr denotes stresses ,e strains and C° a real symmetric positive-definite matrix 
with 

4 =c° j (x,y,z)eL°°(n) 

I lere equation ( I ) is taken as the reference model of the beam. Our task is to ap- 
proximate the 3- spatial-dimensional (3-D) equation (I) by a 1 -spatial-dimensional 
(1-D) beam eq. to arbitrary accuracy of the displacement field. Instead of going 
tin oitgh the term by term scrutinizing as in solid continuum mechanics, we provide 
a unified and systematic approach. This will insight the general pattern and prop- 
erties of the 1-D beam eq. 

I he final goal is the capability of modeling repeated truss structure as a 1-D 
beam. The properties of repeated truss structure, though non-homogeneous (i.e. 

p - p(t,y,z», are periodic along the generic axis x. A homogenization process then 
is needed as will be described in section 8. 
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2 Justification of Taylor’s Expansion 


Before applying the Taylor approximation in the next section, let’s justify its ap- 
plicability to our problem first. Let H = L 2 (Q) be real valued Hilbert space and 


U(t,x,y,z) 
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. U 3 - 
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Let the equation of motion be written as, by equation (1), 

pU + CU = f in ft ( 4 ) 

with well-posed homogeneous boundary conditions specified by forces or displace- 
ments in I\ and 

it,(0, x, y, z) € V, Ui(0,x,y,z)eH whereV = W^ft) 

If all coefficients are in L°°(ft) and the strain energy associated with £ is positive 
definite, then there exists a unique solution 

tii € i 2 ((0,T);V) with u\ and in in L“((0, 7 );7f ) 

Proposition 1 

3 H n X„ -* U strongly in [£ 2 ((0, T); V)] 3 

with X n th- € solution of M n X n 4- A n X n = f n j where 
H n -- H n (y, 2 ) A'n = X n (t, x) € [£ 2 ((ffi T); V )j" 

in [L 2 ((0, T )\ )] 3 = Wi 

H„X n — U in Wi 

Xi, and A' n G [L a ((0,T); tf)] 3 = W 2 , with 
H = L 2 (0, L) V — H l and 


M„(x) 

= J J H*{y,z)p(x,y,z)H n (y,z)dydz 

C„(x) 

= J j T*(y,z)C 0 (x,y,z)T n (y,z)dydz 

/n(U) 

= J J H*(y,z)f(t,x,y,z) d ydz 

[CnfJ ) ( ^]w 2 

= [j 4„A„, $ n ]w 2 
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Proof: 

For each \n{t,x,y, z) 

3 u im — ► Ui strongly in L 2 ((0,T);V) 

u im «; strongly in L 2 ((0, T); ?f) 

Uim —> Mi strongly in L 2 ((0,T)\H) 

where U{ m defined in (0,T) is analytic function in y and z. The strong convergence 
is guaranteed since analytic functions are dense in L 2 ((0,T); V) 

Let 

uM(a,b) = Y, Diu(a)bi 

i=y,z 

« ( 2 \a,b ) = ^ y" Diju(a)bibj (5) 

*=!M 

« ( 3 ) (M) = E E E Dij' k u(a)bibjb k 

*= 2 M j=y,z k= y> z 


By Taylor's theorem 

n j 

«im = 2J jfct «i ( m J ((<,*. 0 1 °),(t,*,y^)) + ^.O.r. (6) 

= J5r„(y 1 «)X im „(* ) *) + /f.O.T. 

By strong convergence of the Taylor series 

H n Ximn — m j m strongly in i 2 ((0,T);V) 

( H n Xim n y — > u- m strongly in L 2 ((0,T);?f) 

HnXimn -* u im strongly in L 2 ((0,T);H) 

Therefore, in general, we have 

H n X n ^U strongly in [L 2 ((0, T)\ V)] 3 (7) 

where 

H n = H n (y,z\ X„ = X n (t,x) and X n G [L 2 ((0, 7’); P)]" 

(H n X n )' - U' in Wi 
H n A'„ — ► U in Wi 
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Also, 


T n t n — c by (H n X n )' — U l 

therefore, 

[pt + CU, $]w, = [/. <t>]w, V * G [i 2 ((0, 7’); V)] 3 

=> -[pf/,$]w, + C*]w, = [/ ^]w, 

or, 

lim { — [pH n Xn , ^n^n]wi + f A' ' f 3»] Wi } — [/, ^n^n]Wi 

ri 

lim {-[//>//„ X„ , * n ]w a + [^f-*o7nf A- - 4]w,} = [H*f, #«]w 3 

n 

Using M n (x), C n (x)andf„(t , x) defined before, we have 

lim { — [Mn X n • 'i'njw’i + \pn*'x ’ e + ]w 2 } = [fn ■ ^n]w 3 
ri 

lim [M„ A‘„ + A„Xn,^n\w, - [/.. . *»]w a (#) 

n 

Liquations (5) and (6) imply 

H n X n — * U with A' n the solution of M n X„ + A ri A u = f n 
Phis completes the proof and justifies the applicability of laylor s expansion. 


3 Taylor Series Approximation 

Let cross sections with concentrated forces be taken as boundary sections (boundary 
points in beam equations) and local effects of applied forces be neglected. Assume 
the physical displacements U, V and W are analytic in y and z so that the 
Taylor series expansion is applicable. We apply the Taylor series expansion, using 
equations (5) and (6), at each cross section x and any time t to have 


U(t.x,y.z) = V + 



8U OHJ f 

dz Oy 2 2 



2 


+ 


d 2 u 

dydz 


yz + H.O.T. 


(9) 


where all terms on the RHS are evaluated at (t t x,0,0). Similar equations can be 

written down for V(t,x,y,z) and W(t,x,y,z). 

The displacement field can well be approximated by a few dominant terms for 

most physical beams. The generalized displacements of the beam eq. can be chosen 
by order of magnitude analysis. For example, for an 8- generalized-displacement 
(8-d) beam eq., we choose 
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We will call X the generalized displacements . 

The approximation of the displacement field up to the specified accuracy then is, 
by equations (9) and (10), 


U(t y x,y,z) = u - <f> 3 y + <j> 2 z + uyz 
V(t,x,y } z) = v + (€23-^1 )* 
W(t y x> y,z) = w + (c 2 ;j + <f>i)y 


( 11 ) 


or 


U 

V 

w 


= HX 


( 12 ) 


with cross sectional shape function H=II(y,z) found from equation (11). The equa- 
tion (11) is known as kinematic assumption in Structural Dynamics, viewed as a 
polynomial approximation to the displacement field. 

Since dynamic eq. of thebeamis completely characterized by mass and constitutive 
properties, the approximation of 3-D eq. by the 1-D beam is equivalent totransforming 
the point properties to sectional properties, ie, from mass density p and constitutive 
matrix C° to mass matrix M and stifTness matrix C , respectively. The mass inertia 
(M) and stiffness (C) have additive property at any given cross-section; therefore we 
can find M and C by approximating the cross-sectional kinetic and potential energy, 
respectively. This additive property justifies the validity of domain extension from 
material-domain to structure-domain. 


4 Energy Equivalence Method 

We can find the sectional mass matrix by approximating the kinetic energy using 
equations (10) and (12). The kinetic energy of a piece of the beam (between any 
two cross sections) is 

KE = ^jj j p(U 2 + V 7 + \V 2 )dAdx 
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i * J J J pX T H T HXdAdx 
- l - J X T MX dx 


It can be shown that, for anti-symmetric mass distribution, 

M = j J pH T H dA - {m,;(x)} 

m X i • • * ~ m b6 


m u 


mn 


mss + 


m 5 5 m 56 m 57 

m 6 6 ^67 

m 7 7 


ni66 “ ^55 


m 55 -f m 66 


( 13 ) 


The M in eq. (13) is the most general pattern of mass matrix for anti-symmetric 
anisotropic 8 -d beams. The first (6 x 6 ) part of M is the general pattern of mass 
matrix for Timoshenko beams. 

We can also find the sectional stiffness matrix C by approximating the potential 
energy of the beam. From equations (2) and (11), we have 


hx 
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(14) 


where T=T(y,z) is a (6 x 9) matrix ; 
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The potential energy of a piece of the beam is 

PE = i j j J e T C°edAdx 

= l J J J c T T T C°Tt dA dx 


\J 


c T C(dx 


The stiffness matrix for anti-symmetric of c 0 - is 


c = J jT T C°TdA = {c tJ (x)} 


Cl3 



C 17 


Cl9" 
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^88 

c 99- 


(16) 


The stiffness matrix C in eq. (16) is the most general pattern for anti-symmetric 
anisotropic 8-d beams. The first (6 x 6) part is that of Timoshenko beams. Knowing 
the general patterns of M and C is very useful, especially in assuming the model 
structure in system identification. Since most of the truss structures ever built are 
at least cross-sectional anti-symmetric, we consider this case only hereafter. 


5 Extension to Lattice Structures 


For non-solid beams, we need to apply the concept of domain extension, from 
material-domain to .structure-domain, so that the results in the above sections can 
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be applied. Let’s take a rectangular lattice structure as an example. The space 
physically occupied by the structure material is called the material-domain. The 
smallest simple-connected rectangular space enclosing the structure is called the 
structure-domain (Q), which includes the space not occupied by the structure 
material^). 

A displacement, field is assumed for the space not occupied by the structure 
material so that the displacement field on structure-domain is in W 1 ^). Therefore, 
the Taylor series expansion and energy equivalence method for calculating sectional 
properties can be applied directly. The sectional properties shall not be affected 
by the introducing of the displacement field in , since both the kinetic energy and 
potential energy arezero in Q e - We can then pretend wearedealing with a solid flexible 
structure in regular shape. 

6 Generalized Beam Equations 

The governing eq. can be found from integrating by parts of potential energy. 

2 PE — j e T C(x)edx = J ( T F dx 
= J(KX' + GX) T F dx 
= J{[-KX} t F' + ( GX) t F } dx + (KX) r F\' 0 
= J X t [—K t F' + G T F } dx + X t (I< t F) |J 

The dynamic eq. in force-acceleration form is 

M(x)X — K t F' + G T F = 0 (17) 


1 T 

Q'n 

Q\ 3 

m(x)x = a' t f'-g t f= m4 M ‘ Qi3 

M' + Q l2 

M' ri - M., 

_ - 2Q 23 

The above two equations are valid for beams which are nonhomogeneous along the 
generic axis. 

From equation (15), 

F = Ce = C(KX' + GX) 
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we have, for a special case of homogeneous beams. 

-I< r F' + G t F = -K t CKX" - (K t CG - G r CK)X' + G r CGX. 
The dynamic eq. in mass-stiffness form is 

MX - K t CKX" - (I\ T CG - G 1 CK)X' + G T CGX = 0 

where 
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with ntij from equation (13) and dj from equation (16) . 

The end boundary conditions are either X{ = 0 or F, = 0. There are 2 
possible combinations, theoretically. For example:*), = 0 for clamped end, K F b - 
( K T CK)X' b + K T CGX b = 0 for free end. 

If there are lumped masses, then the conditions become : 

(a) for interior points 

M b X b - (K t CK) AX' h = 0 

where 

AX b = X( b+)-X(b~) 

(b) for exterior points 

M k X b + sig[{K T CK)X' b + ( K T CG)X b ] = 0 

where 

J 1, for positive ends (ends with positive outward normal); 
st 9 “ \ — 1, for negative ends. 


and 


"■ ' • * c 13 -C12 * 2 c 19 ’ 

• * • • c 23 —C22 • 2c 29 

• ■ ‘ * C33 - C 23 ' 2c 3 9 

'<*cc= ; ; ; : : : :: : <*> 

• C68 

■ ■ • ■ C37 -C 27 C 78 2079 

_ ' C88 - 


7 Timoshenko Beams 

The Timoshenko beam eq. is obtained by deleting the last two generalized dis- 
placements (ie., u and C 23 in X) in equation (19), to have 

M t X t - C t X" - A x X[ + A 0 X t = 0 (25) 

where 

mn ■ 

m 2 2 • 

m33 

m44 

s y m mss mse 

m 5 6 me 6 
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C 1 03x3 

.03x3 C 2 


0 3x3 C\L\ 

03x3 03x3 


A x = L - L t 

’03x3 * ‘1 

Aq = * _ 03x3 03x3 

C 33 C 23 O 3 X 3 C\L\ 

-sym c 22 J 

The force boundary conditions are 

M b X h - C, AX b - L AX b = n (26) 

and the geometric boundary conditions are = specified value. 

8 Multi-scale Averaging Method 

For Periodic Beam-like structure, M and C are periodic in x with period t. 

M = M(x), C = C(x), K and G are constant matrices 

The equation of motion of Timoshenko beam from equation (17) in section 6 can 
be rewritten as 

MX - ( K t CKX ' + K T CGXY + G r CKX' + G t CGX = / (27) 

Let 

B\ = K t CK Bo = K t CG B 3 = G t CG 

We have 

MX - (BiX'Y - ( B 2 xy + B% X' + BsX = f (28) 
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Let 


( 29 ) 


AX = -(BiX'Y - ( B 2 xy + BjX' + B 3 X 
We will consider the following structure 

MX + AX = f (30) 


X(0) = 0 , K t CKX' + I< t CGX = B] X' + B>X = 0 ,at x=L 

.4 = A\ A > 0 

Lets=f£fl, C l {x) = C(f) = C(«), B{(x) = Bi(«), 

B$(x) = B 2 (s), Bi(x) = Bs(s) 

For Xf(t,x) = Xj(t,x, f) = Xi(t,x,s) 


X' = 


dXi , p -i dX < 

Ox ds 


(31) 


We have 


AXi = -(M)' - (B 2 x i y + Bjx' + B 3 X t 
= £~ 2 AoX t + £~ l A\Xi + £° A2X{ 


where 


_ dBy BXi 

A jAj — ■ •: ~~z 2B\ „ ~ 

ds dx ox as 


d . _ dx, 

% 

d 2 x. 


Af)vYi “ ds 1 ' 81 ds 


(B B T ) ® Xi - ^ 8l X 
- ifh - B 2 da X, 


MXi = ~ B] lB~ ~ {Fh ~ B ^Jx + BjXl 


Let [3, 6] 


X l (t, x) = A'o + tX\ 4- P X 2 + PX 3 + ■■■ 

Xi = Xi(t,x,s) i = 0,1,2,--- periodic in s 


(32) 

(33) 

(34) 

(35) 


/ = AX + MX 

= e~ 2 {AoX 0 + (A 0 Xi + PA 0 X 2 + PA„X 3 + • • ■) 
+ r'(yiiXo + (AiXi + PA x X 2 + e 3 A l X 3 + • ••) 
+ {A 2 Xq + IA 2 A'l + Pa 2 X 2 + e 3 A 2 X 3 + ■•■) 

+ M (A'o + iXi + PX 2 + Px 3 ) 

~ C~ 2 (Aq Ao ) 4" ^ 1 (j4i A' o + AoAj) 

+ (A 2 Xo + A\ A i 4- AoX 2 + M Ao) 

+ e{A 2 Xi + A\ A 2 + A 0 A a + M A, ) 

+ <“(•••) 
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The above is valid for all £. Thus, we need 


AoXq — 0 

A x Xo + AqA i — 0 

MX o H- A 2 X 0 + A\X\ + A 0 A 2 = / 
MX\ + A 2 X\ + A\X 2 4 - j4q A 3 = 0 


(36) 

(37) 

(38) 


Proposition 2 
Proposition 3 

with 


A r o = Xo(<,a?) ( not depend on s) 


& A 

A 1 (£, x, s) = Yi y 2 A 0 + Xi(t , x) 


AoYi = 

A 0 Y 2 = 


dBi 

ds 

dlh 

ds 


(39) 


(40) 


(41) 

(42) 


4i and 4 2 are periodic in s since Aj is. Moreover, Yj and Y 2 can be independent 
of t and x and unique up to a constant additive. Note that solution of Vj and Y 2 
are guaranteed since 


l 

L 


dB 

ds 

dB 2 
ds 


~ds = 0 


ds — 0 


Proposition 4 


f( x ) - joj f - B { ~ )ds 
Pi O 151 ds 


d 2 X 0 

dx 2 


\s\ 

I 

\S\ 

W\ l (Bz ~ Bl ^ )dsXo + W\ L M(s)ds Xo 


L f (Pi rT R dYl , nT^l u d. 

PI J\s\ { 2 ~ 02 ~ Bl ~di + 02 Ts“ )ds 1 


dXo 

dx 


A h X 0 + M h X 0 


(43) 


The above is the homogenized eq. found by the multi-scale perturbation method. 
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Proof: 

From eq.(38), to render the sol. of X 2 , we need 

f (^ 2^0 4- A\X\ + MXo)ds = f f(x)ds = \S\f(x) 

J\s\ J\s\ 


( 44 ) 


^ 2 X 0 

A x Xx 


d 2 X 0 T 

~ Bl lW 82 > 

dBj dX A R d 2 Xj 

ds dx 1 dxds 


+ B 3 Xo 

ox 

- (Ba - Bl) 


dXi 

ds 


dBif y d 2 X 0 

' ds 1 1 dx 2 


OX OX 


dY 1 d 2 Xo dY 2 dX 0 , 


ds <9x 2 


ds dx 




dYi 8X 0 dY 2 


ds dx 


ds 


X 0 ) 


ds ax 


dB _ 2 
ds 


(Note that terms associated with and Xi(x) are 0 after integration.) 




+ 


^Ya + 2B,^+(B 2 

ds ds 




dS 2 

ds 


Yi 


ax 0 

dx 


+ ( (« 2 - s ?’) dV ’ 2 


„s + ^ V 0 X » 


jy, ,s’x, <iy, 
(S ' * 


2 ds 7 <9x 

Collecting terms together, we have eq.(43) immediately. 


T dY x dXn RT dY2, x 
Bo + (-£>2 


9 Properties of Homogenized Operator 

A/j, = rii / M(s)ds 
Pi d|S| 

M(s) > 0, M*(s) = M(s) V s 

[M h X, X] = ± [ [M(s)X, X]ds >0 VX ± 0 

Pi J\s\ 

=> M h > 0 


(45) 
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Ml 


= f|l f M‘(s)ds = ^ f M(s)ds = M h 
1*1 J\s\ 1*1 J\s\ 

Mh self — adjoint 

Consider the special case where yY(2,:r) = Xh(t) 

KE C = j [M h X h ,X h ]ds=\S\[M h X h ,X h ] 
•'|S| 

= KE d = [M d X d ,X d ] 


Let 

Then 


X d = T m X h 
Mh ~ ( T m M d T m ) 


(46) 

Therefore, equivalent mass matrix found by Noor’s method agrees with that by 
homogenization theory, though not the case in the stiffness matrix. 


dBi 

ds 


(47) 
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By eq.(41) 
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or 

l 

■ <f c/yj 


=> / [B\~- _ b u ~]ds = 0 

/|S| 
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j [ 5 . 

J\s\ 


ds 


Take # = Y\ 
from eq.(47) and (48) 




(48) 


(49) 


which is positive-definite and self-adjoint. 


B- 


|5| J\s\ 
1*1 J | 5 | 


(50) 
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I (W‘<B, - B.f M.(/- §)] - - £).<* - B, ^)] ) * 
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l.v 


f , <IY; 

U H ' -z - 


1 f t dY-i y , 

lh)' d ?]d* = 0 V'l' € {//' : periodic} 


l akr ^ — Y->, we have 


i=° 


or 


/• , f/Y-> dV 2 

/ Sr'KAi V -* 2 ).*! 

i| 5 | 


ds 


/|5| 
Therefore, 


L , f - *>.<* f - *» = .(.><; + *■ T' ■ W 
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l 5 1 J|5| 

Let L 2 = Z>2(0i M 

[AhXo, A' o]l 2 = j|j J [BT lV ’ V ]i ,, d8 

+ 4i I [(Bz-B^B^B 2 )X 0 ,X 0 }d S 

PI J\S\ 

A h = At Ah > 0 


(52) 


Proposition 5 


/(*) 


-Si 


c> 2 A'o 

<9.c 2 


-S 2 


axo 


-b Xo -+■ MhXo 

= A h A' 0 + M h X o (53) 

where B x S 2 S 3 M h are constant matrices given in eq.(49) (50) (51) (45), respec- 
tively. Also, A h and M h are self-adjoint and positive- definite. 

All the above procedures are formal. In general, B’s are not differentiable and the 
differential eq. should be interpreted in the weak sense. Note that it can be shown 
that the homogenized operators do not depend on B.C. 
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10 Calculation of Operator’s Parameter 

By eq.(4 1 ) 
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II 
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(Vi periodic) => / — -rfs 
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= 0 

Therefore, 



— B\ — C i (Bi >0 C\ a constant, matrix) 


/< 

i|.S| 


Fi periodic => / (7 - )7s = 0 


By eq.(43), therefore, 
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From eq.(42) 
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l |S| 

f Bt l ds)-' 

J\s\ 

nil 

' (Bi — B\ ~p )‘ 

5 | ds 

Ci 

Ci ds 
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<1SU 

f B-'ds )- 1 

|S| 

f,l>0 ’ iu 

= s r 
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•'121 J\S\ 
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J 15 
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b 2 


d r +B ^ d -r )ds 

ds ds 


Using the above results and eq.(43) 

= t4 / (B2 
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^ PI 
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Proposition 6 

Proof: 


PI J\S\ 

B 3 = £3 and # 3 > 0 


C > 0 => [C(KZ\ + GZ2), (K Z\ + GZ 2 ] > 0 ^ Zu %2 
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Take Z\ — — B\ 1 B 2 Z 2 


B\ B 2 
Bj B 3 


Z 1 
Z 2 


Zi 

z , 


= [(B 3 - BjB^ l B 2 )Z2,Z 2 \ 


^ B 3 - Bj Bi 1 B 2 > 0 Vs QE7A 
Proposition 7 77te homogenized eq. is 

e a ! -v 0 6 a-vo 
/w = “ B > - B ' a 7 

+ B3 Xq + MhXo 

= A^Xq T MkXo 


(58) 


where B\ B^ B 3 can be calculated directly from the original B’s parameters as in 
eq.(55), (56) and (57) , respectively. Also , Ah is self-adjoint, and semi-positive 
definite . 

Proof: 

This is a similar result to the Proposition (5). Here we use the numerically calculable 
equations to prove it. The equations for BA have been shown already and the self- 
adjoint and positiveness are as follows: 
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1 

I SI J\s\ 

dX 0 


+ hi / (B 3 -B?Bi'B 2 dsX 0 ,X 0 ] > 0, V An 

Pi JI5I 


(59) 


and 


[A h X 0 M = [- 




— [A'o, —B\ 


d 2 Y 0 

dx 2 

d 2 Y 0 


[Xo,-BI^J + B* 2 ~ + B- 3 Y 0 ] 


x.dYo 

dx 


dx 2 _jB2 ^ +5aV ' ol 


= [A'o, V'o] self-adjoint 


(60) 
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11 Conclusion 


The differences between our approach and the direct averaging method [1, 2] are 
two fold. 

• We know that the stiffness is additive at each cross section and the compliance 
is additive along the generic axis. The energy equivalence method averages 
the stiffness over one repeated cell and thus violates the additivity principle. 
Our approach intrinsically follows the additivity principle. 

• Our approach finds the periodic governing equation first then homogenizes it. 
Namely, we replace the real structure by a 1-D periodic one, then average 
the four matrices to replace it again by a homogeneous beam. The direct 
averaging method averages the properties then finds the governing equation. 
Namely, the method averages two matrices for replacing the real structure 
by a homogeneous beam, then finds its governing equation. 

These are t he major reasons why our approach is more accurate than any previous. 
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1-D HOMOGENEOUS MODEL 


Figure 2: Schematics of Approximation 
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